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Changing  The  Way  Engineering  Is  Conducted 
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Analysis  of  Coaxial  Injectors 

Using  CFD++ 

Presenter: 

Henry  Vu,  Ph.D. 
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FA9300-06-D-0002  and  FA9300-1 0-C-4002 
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Overview 

•  Company  Background 

•  Current  CFD  projects 

•  Coaxial  Jet  Flow  with  Variable  Density 

•  Coaxial  Particle  Laden  Flow 
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Advatech  Pacific,  Inc. 
Background 


•  An  Aerospace  Engineering  Research  &  Development 
Company  Founded  in  1995  primarily  focused  on: 

■  Aerospace  Vehicle  Physics-based  Modeling,  Simulation  and 
Analysis 

■  Electronic  Communications  System  Interoperability 

■  Aerospace  Engineering  Design  and  Analysis  Services 
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Contract  Objectives 


•  Objectives 

■  Modeling  and  analysis  of  fluid  flows  and  heat  transfer  in  support 
of  experiments  performed  at  AFRL  Edwards 

-Complement  experimental  data  by  providing 
detailed  visualization  of  fluid  flow  and  thermal 
distributions  inside  experiments 

-Supplement  experimental  data  by  generating  data 
at  test  conditions  not  performed  in  experiments 

■  Provide  independent  verification  and  validation  for  CFD++  code 
development 
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Contract  Objectives 


•CFD  Projects 

■  Coaxial  Particle  laden  flow  dynamics  to  assist  in  design  of  an 
experimental  apparatus 

■  Coaxial  gas/gas  injector  flow  analysis  for  rocket  fuel  injector 
design 

■  Mixing  of  jet  in  cross-flow  for  pre-burner  studies 

■  Evaluation  of  pipe  flow  conditioning  devices 

■  Conjugate  heat  transfer  studies  in  pipe  flow  for  experimental 
design. 
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COAXIAL  JET  FLOW  WITH 
VARIABLE  DENSITY 

Collaborators:  S.  Alexander  Schumaker,  Ananda  Himansu, 
Stephen  Danczyk,  Malissa  Lightfoot 
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Motivation 

•  Liquid  propellant  rocket  engine  injectors  characterized  by  low 
velocity,  high  density  inner  jet  and  high  velocity,  low  density  annular 
jet. 

•  Conservation  equations  linked  through  density  variations. 

•  Need  to  determine  accuracy  of  RANS  predictions  in  flows  with  large 
density  differences. 

•  Validate  Metacomp  CFD++  tool  for  relevant  flows. 
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Michigan  SEI  Case 


•  PLIF  concentration  data  for  single  element  injector 

•  Low  speed  flow  for  ease  of  modeling 

•  Confined  flow  with  nozzle  at  the  exit 

•  Inner  fluid  (I)  is  air  seeded  with  acetone 

•  Outer  fluid  (e)  is  helium 

9.920'' 

•  Pi/Pe=7.5 


2.000" 


□  .50D" 


Inner  Fluid  (I) 

Air  with  3.35%  Acetone  by 

Volume 

Molecular  Weight  (i) 

29.95 

Kg/kmol 

Molecular  Weight  (e) 

4.003 

Kg/kmol 

Chamber  Pressure 

543829.7 

Pa 

To  (e) 

293 

K 

To  (i) 

293 

k 

po  (e) 

0.894 

Kg/m3 

po  (i) 

6.685 

Kg/m  3 

Uavg  (e)  in  tube 

63.559 

m/s 

Uavg  (i)  in  tube 

12.587 

m/s 

Mdot  (e) 

1.50136E-3 

Kg/s 

Mdot  (i) 

5.93253E-4 

Kg/s 

0.14Si 

0,03^ 

j 

0.059" 


s.ooo^ 


1383^ 
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Michigan  SEi  Case 


Mesh: 

2D  axisymmetric 

93386  quadrilateral  cells,  y+=1 


0.04 


_0  02  Lj ^ i— I— I ^ ^ lJ-j ^ \ lJ-j ^ I lJ-j ^ I lJ ^ ^ i ^ Lj ^ l 

-0.1  -0.05  0  0.05  0.1  0.15 

X 

Case  Conditions: 

•Two  species  (He,  Air/Acetone) 

•Base  Equation  Type:  Compressible  Real 
Gas  Navier-Stokes/Euler 
•Equation  of  State:  Ideal  Gas 
•Turbulence  Simulation:  RANS,  realizable  k 
eps  or  SST 

•Turbulence  Intensity:  2%  or  50% 

•nozzle  geometry  at  the  exit  remove  and 
replaced  with  back  pressure  imposition 
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Michigan  SEI  Case 


Centerline  Mass  Fraction 


Measurement  k-eps  k-eps,  50%  inlet  turb  intensity  SST,  50%  inlet  turb  intensity 


•At  least  15%  error  in 
concentration  along 
the  centerline  in  the 
near-field 

•No  improvement  by 
increase  turbulence  at 
inlet 

•SST  turbulence 
model  makes  result 
worse 


Is  error  due  to  problem  in  modeling  scalar  or  velocity  field? 
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Validation  Case 


Ld 

CH 


•2D  axisymmetric  compressible  RANS 
•Density-based  solver 
•Realizable  k-eps  turbulence  solved  to 
the  wall 

•T=300  Kand  P=101325  Pa 
•Center  inlet  boundary: 

•He  normal  velocity  =  24.45  m/s 
•Air  normal  velocity  =  10.5  m/s 
•CO2  normal  velocity  =  9.0  m/s 
•Outer  inlet  boundary:  Air  with  normal 
velocity  of  0.9  m/s 

•Outlet:  Simple  back  pressure  outlet  of  0 


Pa. 


Gas 

Uj  fm/sl  (Expt) 

Uj  [m/sl  (CFD) 

Re, 

p/p. 

He 

32 

31.4 

7000 

0.14 

Air 

12 

12.68 

21000 

1 

CO, 

10 

10.75 

32000 

1.4 

Amielh,  M.,  Djeridane,  T.,  Anselmet,  F.,  &  Fulachier,  L.  (1996).  Velocity  near-field  of  variable  density  turbulent  jets.  Int.  J.  Heat  Mass  Transfer , 
2149-2164. 

Djeridane,  T.,  Amielh,  M.,  Anselmet,  F.,  &  Fulachier,  L.  (1996).  Velocity  turbulence  properties  in  the  near-field  region  of  axisymmetric  variable 
density  jets.  Phys.  Fluids  ,  1614-1630. 
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Validation  Case 


•234880  quadrilateral  cells  with  a  y+  of  1  or 
less  at  all  wall  locations  except  the  outermost 
wall. 

•Run  up  length  =  1.5  m 
•Total  domain  length  =  2.7  m 
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Validation  Case  Results 


He 


CO 


2 
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Validation  Case  Resuits 


1.2 


O 

"^o 
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Radial  Profiles  of  Mean  Mass  Fraction 
- 1 - ^ ^ ^ - 1” 
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1.2  I - r 


O 
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Radial  Profiles  of  Mean  Mass  Fraction 
n - 1 - ^ ^ 


1 - ^ 
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0.2  i - ^ 
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X/D,=0.20  CFD 


□  X/Dj=0.20  Expt 
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Validation  Case  Results 
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Validation  Case  Results 
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Validation  Case  Results 


Air 
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Validation  Case  Results 


1.2 

Radial  Profiles  of  Mean  Longitudinal  Velocity 
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Validation  Case  Results 


Radial  Profiles  of  TKE 


X/Dj=0.20  CFD 
X/Dj=0.20  Expt 
X/Dj=5.00  CFD 
X/Dj=5.00  Expt 
X/Dj=10.0  CFD 
X/Dj=10.0  Expt 
X/D=15.0  CFD 
X/Dj=15.0  Expt 
X/Dj=20.0  CFD 
X/Dj=20.0  Expt 


Radial  Profiles  of  TKE 


X/Dj=0.20  CFD 
X/Dj=0.20  Expt 
X/Dj=5.00  CFD 
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He 


Radial  Profiles  of  TKE 

0.4 
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1  1  1  1 

1  1  1 
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1  1  1  1  1 

0.35 
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Validation  Case  Results 


He 


Velocity  Half-Width 
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Validation  Case  Summary 


•  Case  exhibits  same  deficiency  in  centerline  mixing  for  He/Air 
in  the  near  fieid  whiie  improving  in  the  far-field. 

•  Radial  profiles  of  mean  quantities  generally  good,  except  at 
centerline. 

•  C02/Air  case  shows  accurate  mixing  in  the  near-field  and 
worsens  downstream. 

•  Both  CFD  and  experiment  show  low  sensitivity  of  density  ratio 
on  jet  spreading.  Mixing  occurring  by  entrainment  of 
surrounding  fluid. 

•  Turbulence  is  generally  over  predicted  everywhere  except  in 
the  near  field  for  He/Air. 

•  Deficiencies  in  centerline  mixing  may  be  correlated  to 
turbulence. 
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Future  Work 


•  Try  difference  turbulence  models. 

•  Modify  k-eps  to  compensate  for  low  near-field 
turbulence. 

•  Try  Large  Eddy  Simulation  (LES). 

•  Proceed  to  more  relevant  flow  conditions  for  rocket 
injectors. 
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COAXIAL  PARTICLE  LADEN 
FLOW 

Collaborators:  Ananda  Himansu,  Alireza  Badakhshan,  Stephen 
Danczyk 
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Motivation 


•  Experimental  set  up  in  lab  to  evaluate  novel  fuel  ignition  strategies 

•  USI  device  with  a  shroud  to  induce  swirl  in  coflow 

•  Use  CFD  to  help  design  apparatus  and  flow  conditions 

•  Some  earlier  CFD  work  helped  in  redesign  of  shroud  for  better 
swirling 

•  A  parametric  study  to  test  the  effects  of  swirling  flow  showed  some 
unusual  results 

•  Need  to  validate  dispersed  phase  modeling  features  in  CFD++  to 
confirm  results  and  determine  important  conditions  to  accurate 
modeling 
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Eulerian  Dispersed  Phase 
Modeling  in  CFD++ 


Dispersed  Phase 

Mass: 

^  +  V  ■  (ppiUpi)  =  mi 

Momentum: 

^  ^  ^  +  V  ■  (ppiUpiUpi)  =  FDi  +  Fv.M.i  +  Ft.d .i  -h  FLi  -h  Fp,G.i  +  Fsi  +  miUpi 

Energy: 

^  ~  Conduction  source  term  +  Radiation  source  term  F^i^pi 

—  ^^pi  F  "^i^pi 

Number  Density: 

8N- 

^  +  V.(JVifl,0=O 

Melting  Fraction: 

^  ~  Melting  source  term 

Material  ^  ^  t  ~  /  c\ 

Density:  ~  "I"  Pp^^^nd 


Particle  n  =  (  ^ 

radius:  ' 
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Eulerian  Dispersed  Phase 
Modeling  in  CFD++ 


Continuous  Phase 


Momentum 

source 

terms: 


^  (^^Di  -I-  ^V.M.i  -I-  , 


Energy 

source 

terms: 


NP 

E 

i  =  l 


Ppi  ^pi  /aTi 


(Tf  —  Tpi)  -\-Uf  ■  (F^.  +  Fv.M.i  +  ^T.D.i) 


Relevant  Source  Terms 


Interphase  Drag:  Fd; 

Ppi 


fDz.^  ^  . 
-—[Uf  -Upi)  , 


Turbulent 

Dispersion: 


—r 

i'T.D.i 

Ppi 


/  Vrjpi  _  VVf\ 
Tuz  Plf  \  Vpi  Vf  ) 
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Model  Validation  Case 


Jet  airflow  ^  Jet  airflow 


Liquid  flow 


Nozzle-tube 


Mesh  screen 


Sonotek  nebulizer 


Poly-dispersed 
turpentine  droplets 

d  =  1-90  urn 

=  2.4  m/s 

Low  turbulence 
intensity  of  1 .4% 


Nijdam,  J.,  Langrish,  T.,  &  Fletcher,  D.  (2008).  Assessment  of  an  Eulerian  CFD  model  for  prediction  of  dilute  droplet 
dispersion  in  a  turbulent  jet.  Appl.  Math.  Model.,  2686-2705. 
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Model  Validation  Case 


Droplet  Inlet  Boundary  Condition  (1D  downstream) 


Variable 

Constants 

Excess  axial  mean  velocity  (m/s) 
t/e  =  t/eoLl  -  [sm{{R/Rl/2u)")r\ 

( 1 )  Peak  exeess  axial  mean  velocity,  Ueo 

Uco  =  -0.05148  (1  +  21.97494 

where  d  is  droplet  diameter  (pm) 

(2)  =  3.9320, 112  -  1.222,  4.978 

Radial  mean  velocity  (m/s) 

r  =  0 

Volume  fraction 
r  =  /  q  exp[-.4(y?//?i/2VFrj 

( 1)  Peak  volume  fraction  Iq 

4.32H-08  (5  pm),  6.91E-07  ( 15  pm), 

5.07H  06  (25  pm),  9.50H  06  (35  pm), 

1.47H-05  (45  pm),  2.16E-05  (55  pm), 

2.04E-05  (65  pm),  7.43E-06  (75  pm), 

2.50E-06(85  pm) 

{2)  A  0.6942,  n  2.1543,  7?i/2vf  =  3.9480 

Gas  turbulent  kinetic  energy  (m^s^) 

A  -  Dexp[-.>t(/? -f  C 

(1)  A  =  1.011,  5  =  4.904,  C-  1.75,  Z)=8.73,/?-  1.418 

Gas  turbulent  energy  dissipation  (m^s^) 

~  0.2D 

( 1 )  Z)  0.0098  m 

Nijdam,  J.,  Langrish,  T.,  &  Fletcher,  D.  (2008).  Assessment  of  an  Eulerian  CFD  model  for  prediction  of  dilute  droplet 
dispersion  in  a  turbulent  jet.  Appl.  Math.  Model.,  2686-2705. 
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Model  Validation  Case 


Case  Conditions:  *79000  quadrilateral  cells 

•Two  species  (Air,  Turpentine)  .Total  domain  size  =  0.05  m  x  0.4  m 

•Base  Equation  Type:  Compressible  Real  Gas 

Navier-Stokes/Euler 

•Equation  of  State:  Ideal  Gas 

•Turbulence  Simulation:  RANS,  realizable  k-eps 

or  SST 

•Turbulence  Intensity:  1.4% 

•10  urn  droplet  size  bins  from  5-85  urn  diameter 
•Temperature-based  inlet  profiles  generated  with 
Matlab  script 
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Results 
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Results 
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Results 


R  [mm] 
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Summary 


•  Excess  axial  velocities  of  smaller  droplets  up  to  50%  slower 
than  experimental  while  larger  droplets  show  better 
agreement 

•  Error  for  peak  volumetric  fluxes  higher  for  larger  droplets- 
may  be  due  to  the  fact  that  they  carry  more  volume 

•  CFD  agreement  with  radial  velocities  is  poor-  may  be 
measurement  error 

•  EDP  features  available  In  CFD++  may  be  missing  important 
physics 
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Future  Work 

Parametric  study  of  effects  of  co¬ 
flow  velocity  components  on  droplet 
mixing 
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Future  Work 


Use  CFD  to  design  and  determine  conditions  for 
experimental  setup  for  testing  combustion  ignition 
strategies 
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Questions? 
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